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AB^RACT 

Underwater sound was scattered from a wind-agitated random 
rough model sea surface c The surface energy spectral density, 
the surface rms height, and the surface correlation function were 
measured* The mean square magnitude and probability distribution 
of the specularly scattered sound were measured for various angles 
of incidence between 35° and grazing. The results are compared 
with the basic theories expressed by Eckart and Beckmann and 
shadowing corrections derived by Wagner. The frequency spectrum 
of the scattered sound was studied and found to contain doppler 
shifted components, contrary to theory which predicts no doppler 
shifting in the purely specular direction. A detailed analysis 
of the origin of this doppler shifting is presented, and it is 
shown that interaction between axial and off-axis radiation is 
the explanation of the observed frequency broadening. 
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Specular Scatter of Underwater Sound 
fron a Wind-Driven Model Sea Surface 
1 , Intr oduct ion ; 

Backgr ound 

When underwater sound is incident on the surface of the water, 
it is scattered* This scattering is dependent upon the character- 
istics of the surface at the time* If the surface is smooth, the 
scattering will be pure specular reflection; that is, angle of 
incidence equals angle of reflection; both are measured with re- 
spect to the normal to the smooth surface. If the surface is rough 
the scattering is no longer purely specular, but energy is scat- 
tered in all directions* Two directions of scattering are parti- 
cularly interesting, back scatter and specular scatter* 

In recent years, since the mid 1950's, a great amount of work 
has been done in studying the back scattered energy. This case is 
relatively easy to study at sea, and many experiments have been 
performed. In 196? Ball and Carlson studied this problem in the 
laboratory at the Naval Postgraduate School under controlled cir- 
cumstances. They used a water surface which was roughened by 
blowing air over it. Using this reproducible surface they examined 
the relative mean acoustic pressure and doppler shift of the back 
scattered energy. Their experimental apparatus will be referenced 
throughout this report. 

Distribution of the Scattered Sound Pressure 
Our experiments were conducted to measure the effects of water 
waves on scattering of underwater sound in the specular direction. 
The applicable theories on the subject are based on the Kirchhoff 
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Assumption, which asserts that plane reflection coefficients can be 
used on a point by point basis for a gently curved surface. The 
basic theories assume that there is no shadowing or multiple scat- 
tering. Wagner (I966) has calculated a correction for the shadowing 
assuming ray propagation. This correction becomes important near 
grazing incidence. Our model experiment includes an oceanographic 
study of the scattering surface and permits a check of both the 
Kir chhof f-based theories and Wagner ^s correction. 

The experiments to measure the distribution of the scattered 
sound pressure were conducted by LCDR scheible and the results were 
presented at the May 1968 meeting of the Acoustical Society of 
America in Ottawa. 

Spectrum of the Scattered Sound 

Observations of scattered sound pressure have shown it to be 
both amplitude and frequency modulated. The most recent attempts 
to describe this behavior have been made by L. Wetzel (1966) and 
B.E. Parkins (1967)# Although Wetzel ^s descriptions are explicitly 
derived for long EM waves, his conclusions are directly applicable 
to the acoustic problem. Wetzel bases his developments on a re- 
sonant interaction between the incident energy and an appropriate 
frequency component in the surface. Parkins uses physical optics 
approximations to describe the reflections occur ing at the surface. 
Both of these developments arrive at predictions of doppler shifting 
of energy upon scattering ftom the surface. They predict doppler 
shifting for all directions of scatter except the specular. They 
thus imply that the energy scattered in the specular direction would 
be constant in amplitude and that it would not contain any frequencies 
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other than the incident frequency. However, it has been observed 
that specularly scattered energy does fade and may contain frequency 
or phase modulation. Mellen (1964) has observed this frequency be- 
havior for back scattering. The frequency spectra of specular 
scattered sound Jias been measured at sea, near Bermuda, by Roderick 
and Cron (1968). * Their observations show the frequency spread to 
be less than 1 Hz under all ocean conditions and less than .2 Hz 
under normal conditions. 

The experiments to measure the spectrum of the scattered sound 
pressure were conducted by LT Fowler and the results were presented 
at the May 1968 meeting of the Acoustical Society of America 
(Supplementary paper D3A) . 
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2 . Theory : 



Spectrum of the Scattered Sound 

The following treatment is a summary of Parkins (1967)* **When 
a monochromatic acoustic wave is incident upon the time varying 
surface of the ocean, the envelope of the reradiation varies in 
time or fades.” Parkins claims that this fading is a result of 
doppler shifting of the scattered pressure. His contention is 
that the incident radiation is ”locally shifted in frequency” by 
scatter ers moving on the surface. He assumes the surface to be a 
stationary, homogeneous, Gaussian process that is described by a 
thr ee-dimens ional spectrum . 

The locally incident pressure is given by 





•r 
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e 



1 A 



A 'A 

+ jm + kn) 



where 1 , m, and n are direction cosines. Using developments 
similar to Beckmann (1963), the scattered pressure in the far field 
as a function of position and time is given by 
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X = acoustic wave length 




^ dx dy 




area of integation = a square 21 on a side 
ki, k^, and r are defined on Figure 2 . 1 . 
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Fifure 2*1 7>jo Dimensional Ceometry of Scattering From the Surface 
( y direction is out of the paper) 



To find the spectral density of the scattered pressure, Parkins 
transforms its auto-correlation function 

S(Q) = j R(T)e^^^dT 

-OS 

* 

where R(t) = \0^(P>1)0 + t )) , that is the ensemble average, 

which is the time average if you assume ergodicity. He then shpws 

2 2 

that for a slightly rough surface, a « 1 

S(n)= B exp(-K^ a^)^6 (n-ou)l6 
+ a^G(K^,Ky,n^)j 

where B = (0'^/2nR (K-k,/K 

o 1 z 

CD ^ 

G(K^,K ,n-(J)) = J G»(K^,K ,T)e^( dT 

^ -CX> ^ 

CD 

= /y p exp [-i(K^u + k^v)l du dv 

-OS 

U = X^-X 

V = y»-y 

p = surface auto-correlation function as a function of space 
and time 

a = surface rms height 

Assuming a Neumann -Pier s on directional wave spectrum and 
integrating, he then shows the spectral density of the scattered 
sound pressure as 
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S(fi) 



B exp(-K^a^)[ 6(0-uJ) l6 Sin\^^ Sin\^Ji/K\^ 

+ 2X\^a^ [[A»(K^,K^)1^6[n-U)-(K^+Ky)\^] 

+[ A • ( -K^ , -Ky ) 1 ^ 6 [ n + (K^+K^ ) ] 



The importance of this equation is that is shows the reradiated 
pressure to be composed of a specularly reflected component at 
the same frequency as the incident pressure and two components 
that are doppler shifted by 



(k" - 

X 



y' 



Parkins says, "The scattered components are perhaps sur- 
prisingly, given by delta functions*" He further states, "It 
seems likely that the spectrum of the scattered radiation will 
be narrow but not have delta-function behavior. The spectrum of 

the specularly reflected component depends only upon the gravity 

2 2 

wave spectrum through the phase modulation term exp(-K^ o ), and 

it will consist of a single line at Q = (JU as long as the sea, as 

seen from the angles of incidence and observation, is only slightly 

2 2 

rough." By slightly rough he means the roughness factor K a = g 

z 

of Beckmann is less than JL . 

It should be noted that Parkins obtains both positive and 

negative doppler shifts by considering scattering from opposite 

sides of the origin, that is + K , + K . Mellen (1964) measured 

X y 

the spectrum in the back scatter case and found doppler shifting 
in only one direction- Further, his doppler shift is about one- 
half that predicted by Parkins. 
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In 19^6, Wetzel, by considering a resonant interaction between 
incident HF radiation and surface wave lengths, also predicted 
doppler shifts. His scattering geometry is shown in Figure 2.2, 

He developes the relationships 

K_ = + k (Sin 0. - sin 0 Cos 0 ) 

1 — ^ 1 s ^s' 

k (Sin 0 Sin 0 ) 

2 — ^ 3 s ^ 



where = the x component of the wave propagation number 
= 'the y component of the wave propagation numb^er 
K = the wave propagation number 
k = 2TT/X 

X = wave length of the radiation 

He argues that for conical scatter, ^ ^ ^ doppler 

shift is given by 

= (gk Sin 0^)"^ [2(1 - Cos 0g)l^ 

where g = the acceleration of gravity. This equation also predicts 
zero doppler shift in the specular direction. 

Both of these developments have been for plane waves, that is, 
beam pattern divergence has not been considered. The assumption 
of a plane wave front is justified when the surface is far from the 
source and only a small portion of the wave front is considered. 
However, it seems likely that this will not always be true when 
considering sonar energy interacting with the surface. 
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:'igure 2.2 V.'etzel's Scattering Coordinate System 
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Distribution of Scattered Sound Pressure at 



Near Grazing Incidence 

The magnitude of the specularly scattered sound pressure is 
a function of the relative roughness of the surface • The most 
commonly used measure of roughness was defined by Beckmann (1963) 
as the ^roughness factor^^ g, which is defined by the equation 
g^ = kcj (Cos Q + Cos where k = 2n/X is the propagation con- 

JL w 

stant of the sound, a is the rms wave height, and and 0 ^ are 

the angles of incidence and reflect iop« For specular scattering 
2 2 

Q = Q . g = K o in Parkings notation* 

X. ^ z 

An equation for the relative mean square pressure of specularly 
scattered sound, measured relative to mirror reflection, has been 
developed by Medwin (196?) : 






-9 , ^ V 

e + 



A Cos 0 



64 






s(g) 



In that paper S(g) was defined, and it was shown to be approximately 

equal to g for g < 0 * 3 * For our experiment, the second term, which 

represents incoherent scatter, is always less than 2 % of the first. 

-g 

Therefore it has been neglected, and we have used e as the pre- 
dicted value of relative mean square pressure. 

Wagner (I 966 ) derived a shadowing function making use of ray 
geometry and statistical mathematics. For specular scattering this 
funct ion 



S(Q) = (1 - erf v/4b 



where B=^/(3'nv) ^e ^ +(ttv) 



- v"^ 



- erfc V 



j 
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and V = Cot 0/2(a/L). 



a = rms wave height 
L = surface correlation distance 
0 = angle of incidence or reflection 
Wagner also shows in that paper that the correct application of 
this correction is to multiply the predicted value of relative 
mean square pressure by the square of this shadowing function, 

Medwin following Beckmann also developed an equation for the 
variance of the magnitude of the relative pressure for specular 
scatter ing: 



From this equation it can be seen that the variance should depend 
on the sixth power of the frequency and on the fourth power of 
the cosine of the angle of incidence. 




Substituting the definition of g above and c/f for X gives 
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3 • The Model Sea Surface; 



General 

The physical arrangement of the experiment is shown in Figure 
3.1* A large anechoic water tank was covered with plywood 3V* 
above the water surface. Air, from three centrifugal fans, was 
forced between this covering and the water to generate a model sea 
surface. At the opposite end of the tank we have an "anechoic*^ 
beach made out of six layers of small-mesh screen. 

In order to compare the behavior of sound pressure that is 
scattered from the surface with theory, one must know the surface 
energy spectrum, the surface correlation function R(d), the sur- 
face rms height a, and the surface rms slope Tj. 

These developments usually assume a NeumannrPier son surface 
spectrum and a correlation function of either e or e , 

where the correlation distance L is defined as the distance at 
which the correlation is reduced to one e-th of its value at zero 
separation . 

For the surface used in this research the surface energy 
spectral/ dens ity , the spatial correlation function, and the surface 
rms height were measured. The frequency spectrum and rms height 
were measured at several positions on the surface (Figure 3.2). 

The rms slope was previously measured by Ball and Carlson. The 
rms height and slope were calculated from the energy spectral 
density and compared with measured values. 

The experiments that were conducted to measure the surface 
frequency spectrum, rms height, and correlation function were per- 
formed jointly by LCDR Scheible and LT Fowler. 



I 
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Figure 3.1 Side Vievf of Anechoi^ lank 
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I'igure 3*2 Top Viev of Ar»erhoic Tank 



Energy Spectral Density 



Electr onics , Figure 3*3 shows the circuit used to measure 
the energy spectral density E(f) of the surface. The oscillator 
has a stable 10 kHz, 1.5 volt output. A resistive probe (Figure 

3 . 5 ) is placed in series with the output impedance of the oscil- 

; 

la tor. This probe was designed and analyzed by Ball and Carlson. 

i 

t 

As the wetted length of the probe changes, its resistance varies 
and produces an amplitude modulation of the oscillator output. 

This modulation is then an electrical analog of the water surface. 
It is necessary to use at least a 1 volt oscillator output to 
ensure the modulation will be in the linear portion of the de- 
tector (Figure 3*4). Keller ( 1968 ) also used this detector and 
found it to be linear when the input is in excess of .5 volts. 

The detector output is the electrical analog of the surface. 
This analog has about a 1 volt DC component and a 10 kHz ripple. 
The differential amplifier (linear down to DC) is used to remove 

i 

the DC and ripple. The DC is removed by subtracting an equal 
value, which is set on the variable DC supply. The differential 
amplifier has an internal, variable, low-pass filter which is set 
at 200 Hz. The filter is flat from DC and down 3 db at 200 Hz. 
Because the surface frequencies are less than 100 Hz, they are 
passed by the filter but the 10 kHz ripple is blocked. The dif- 
ferential amplifier output is then recorded at .375 ips . The tape 
recorder is used in the FM mode to be sure the frequencies from 
IX to 100 Hz are preserved. The response of the recorder is DC 
to 100 Hz + 1 db, S/N 35 db. The tape recordet is increased in 
speed by 10 , or 100 for play back to increase the surface 
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Figure 3,3 Surface V.'ave Analysis 
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Figure 3.i4 



Detector 










Figure 3.5 Surface Probe for Measuring Wave Heights 
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frequencies to values that can be analyzed. This speed-up is 
necessary because the surface spectrum is predominantly from 0 to 
20 Hz, and the wave analyzer can not accurately measure frequencies 
below 10 Hz. 

The wave analyzer is used to measure the spectrum of the tape 
recorder output. When the tape recorder is played at 3*75 ips , 
the 10 Hz filter is used. In the first experiments performed, a 
complete spectrum was obtained for both replay speeds. It was 
found that the predominant surface energy was at 3.8 Hz and that 
one need only analyze with a replay speed of 3*75 ips. This allows 
accurate spectral analysis of surface frequencies in excess of 1 Hz. 

The output from the wave analyzer is a time averaged voltage. 
With the analyzer set to a given frequency, it was found that the 
energy in that component varied with time. This variation was 
partially removed by using a slow response on the analyzer and 
was eliminated by eye averaging. The analyzer voltage was con- 
verted to a 1 Hz bandwidth by multiplying by the appropriate fac- 
tor (.65 for 3 Hz BW, .375 for 10 Hz BW, and .159 for 50 Hz BW) . 

2 

The analyzer voltage is converted to cm by multiplying by the 

system calibration. The square of this result is the energy den- 

2 

sity at that frequency with a 1 Hz BW (cm -sec). A third correction 
factor is necessary to account for the multiplication of frequency. 
When the surface analog is replayed at 3*75 ips and the wave an- 
alyzer is set for a 3 Hz BW, the real frequency BW is in fact 

.3 Hz. Therefore the wave analyzer output energy, (voltage X BW 

2 

correction X conversion factor) , must be multiplied by 10. Like- 
wise, when the replay speed is 37.5 ips, the energy must be mul- 
tiplied by 100 . 
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Calibration > To calibrate the output voltage as a function 



of wave height, the probe depth was changed in increments and the 
DC output from the tape recorder noted. To verify the linearity 
of the tape recorder, the DC voltage input to the recorder was 
also recorded. Figure 3.6 shows a sample of these results. A 
calibration was obtained for each spectrum obtained. The conver - 
sion factor is equal to 

A voltage DC 1 inches _ A volts rms 

A probe height (inches) 2 . 54 cm cm 

Although the voltages recorded were DC, the displacements were 
also DC; therefore, the results can be expressed as volts rms to 
cm rms. Ball and Carlson verified a static calibration to be 
valid for frequencies up to at least 100 Hz. 

Results . The surface energy spectral density was ^determined 
for stations 1 C, 1 R, and 2 C of Figure 3.2. These spectra are 
presented in Figures 3.7 through 3.9 t Spectra were determined with 
and without the beach installed. It was found that the spectra 
without the beach were higher and peaked to the left of the ones 
with the beach. This is the result of reflection from the end of 
the tank. Since the objective of this report is to compare the 
scattered sound pressure with a sea-like surface, it was desirable 
to have the surface propagation vectors in the downwind direction 
only; therefore, it was decided to use the beach. 

Initially the spectrum was obtained with the probe attached 
to the plywood tank covering. It was felt that this plywood might 
vibrate and introduce noise into the system. Although this was 
not found to be true, the probe mounting was shifted to the metal 
grating which covers the tank. 
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Figure 3,6 Surface Probe Calibration 
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Figure 3»8 Surface Spectral Density ( station 1C) 
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Figure 3*9 Surface Spectral Density ( station 2C 



The system noise level was found to be down 10 db at 1 Hz, 

40 db at 4 Hz and 20 db at 10 Hz. The spectrum of the noise was 
determined using the same procedure as the surface spectrum except 
that the fans were off. 

Examination of the spectra in Figures 3*7 through 3*9 yields 
the following. The energy at station 2 is higher than at station 
1 by about 1 db, which is within experimental accuracy. The energy 
1.5 feet off the center of the tank is down by about 5 db, which 
indicates that the air from the fans was funneled to the center of 
the channel. The effect of the beach was to remove the energy 
reflecting from the end of the tank, thus lowering the total energy. 
All spectra with the beach installed peak at 3-8 Hz, have an f 
dependence near the peak reducing to f * dependence with in- 
creased frequency. The f dependence is in good agreement with 

the Neumann-Pier son spectrum and Phillips ( I 966 ) theory which 
predicts f 

Using the surface energy spectral density it is possible to 
find the rms height a and the rms slope S. These two parameters 
are used in various theories to predict the behavior of the scat- 
tered sound pressure. 



Integration of Figure 3*8 yields ^ o of 0.22 cm, which is in good 
agreement with that found using a sampling technique for the same 
surface as shown later in this section. O.M. Phillips (I 966 ) 
shows that 00 
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ich is equal to (2n)^ J' — L - ) . df , where A is the surface wave 



A 



the values of S ^ , and = S ^ + S ^ 

ud cross ud cross 



length. Integrating Figure 3*10 yields E = .0l8. Using the 

2 2 2 

approximation given in Beckmann, S = 2 a /L , where a is obtained 
from integration and L is the correlation distance (3*2 cm for 
up-down wind and 10,2 cm cross-wind as shown later in this section), 

were calcu- 
lated and compared with the optically measured values of Ball and 
Carlson in Table 3.1. The relatively weaker cross-slopes in the 
present research is attributed to the completely covered water 
surface which maintained the wave build-up predominately in the 
windward direction. The reduced total slope is probably due to 
increased drag due to that greater covering. 



Table 3.1 



rms slopes of the surface 



(given by S^= 2CJ^/L^) 



Scheible 
and 
Fowl er 

Ball and 
Carlson 



S ^ S ^ 



ud 


cr OSS 




.0096 


.0008 


.0104 


.0097 


.0068 


.01655 
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Figure 3»10 Energy Spectral Density Divided by the (Surface Wavelength) 



Analysis of the Energy Spectral Density 
It is believed that the surface height can be described as 
h(t) = E(t) Cos 0(t). If this function could be determined, many 
powerful mathematical tools could be applied to thoroughly analyze 
the surface. That the surface can be described by h(t) is hinted 
by the fact that the energy in a given frequency band is time 
varying. Attempts were made to find E(t) and 0(t) by studying the 
amplitude variation of the energy centered about a given frequency 
f^ and by studying the frequency modulation of the surface. 

Figure 3-H shows the surface analog as measured at station 1 C, 

It is obvious that there is some amplitude and phase/frequency 
var iat ion • 

Amplitude Modulation ^ It is thought that the amplitude 

variations of a given frequency component of the surface 

may have a direct bearing on the scattered sound pressure. To 

find a description of these variations. Various experimental 

attacks were tried. First the wave anlyzer was set to f and 

o 

the filtered output was analyzed. The three available filter 

settings (f^ = 3, 10, 50 Hz) were used. Since the amplitude 

variations were at a very low frequency, the surface analog 

was replayed at 37.5 ips (f^ x 100), Figure 3.12 through 3.l6 

show the time variations of given f ^s and Af ’s; in these 

o o 

figures all times and frequencies have been reduced to real values. 
The output from the wave analyzer was processed as shown in 
Figure 3.17. The output from the second tape recorder was cali- 
brated to the output of the first tape recorder (surface analog). 
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Figure 3*11 Surface Analog ^ ,2 sec/cir.) 




Figure 3,12 Wave Analyzer Output ( Fq = 3*8 Hz,AfQ=.03Hz) 
( ,001 sec/cm) 
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Figure 3.13 VJave Analyzer Output (fQ=3.8Hz,AfQ=,lHz 
.001 sec/cm) 




Figure 3.1U Wave Analyzer Output 

(f0=3.8Hz^f0=.5Hz, .001 sec/cm) 
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Figure 3.15 V«ave Analy'^er Output 

(f(-=ii.8Hz, fQ = .0i2, ,001 sec/cm) 




Figure 3.16 Wave Analyzer Output 

(fQ=2.8Hz, fo=.5Hz, .001 sec/cra) 
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^lfl= 3Hz 



Figure 3.17 Circuit to Measure the Energy Spectral Density 
of a Filtered Surface Component 




Figure 3*18 Circuit to Measure the Wave Analyzer Filter 
Response 
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• 003 Hz in real fre- 



The secona wave analyzer was used with Af^ = 

quency. The calibration was that the energy spectral density out 

2 -5 2 

of the second wave analyzer was equal to (mv) x 1.07 x 10 cm -sec 
referenced to 1 Hz Bw. 

The theory behind this analysis is as follows. If a pure 
tone is amplitude modulated, the modulation spectrum will appear 
symetrically on either side of the tone in the frequency domain. 
Therefore, if one filters f^ and measures the spectrum in the 
vicinity of f^, the amplitude modulation will appear as side bands. 
This is exactly what was done; the results appear in Figures 3.19 
through 3»22. 

To verify that these results were not just functions of the 
filter in use, filter response was measured using the circuit in 
Figure 3.l8. The SG-85 generates a carrier frequency f^ at 10 kHz, 
210 mv. This carrier is externally modulated by the HP 3300A. 

The output from the SG-85 is an amplitude modulated signal. The 
modulation frequency was varied in steps and the output from 

the SG -85 was fed to the first wave analyzer, set to 10 kHz, 

Af^ = 50 Hz. The second wave analyzer was varied to measure the 
output of the first. By holding the amplitudes of f^ and f^ con- 
stant, it was possible to determine the filter characteristics of 
the first analyzer. The result is shown in Figures 3*23 and super- 
imposed on the Figures 3*19 through 3*22. 

Although the theory is correct, it was not possible to find 
the spectrum of the f^ component, because f^ is not a pure tone 
but a component in a band of frequencies. It is not possible for 



44 




45 



Figure 3.19 Energy Spectral Density of the Surface 2,8 Hz Component, .3 Hz Bandwidth 
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Figure 3.20 Energy Spectral Density of tne Surface 3-B Component. .5 Hz Bandwidth 
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Frequency (Hz) 

Figure 3.21 Energy “Spectral Density of the Surface J,8 Hz Component, .L Hz Bandwidth 
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Figure 3*22 Energy Spectral Density of the Surface U.8 Hz Component, .5 Hz Bandwidth 
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Figxire 3*23 Response of the GR 1900A Wave Analyzer Filter ( f = 50 Hz) 



the filter to pass f^ and its amplitude modulation alone; it must 
pass the components of the surface near f^. This is shown by the 
fact that the results in Figures 3.19 through 3.22 are nearly re- 
productions of the energy spectral density near f^. 

The frequencies near f^ that the filter passes tend to mask 
the amplitude modulation of f^. As the bandwidth is increased, 
the amplitude modulation appears to increase (Figures 3.20 and 3.21), 
a result of the added frequencies being passed. To gain additional 
insight into the time-varying surface characteristics a digital 
analysis technique, with its inherent flexibility of filter speci- 
fications, was applied. 

To find the description of the amplitude variations of speci- 
fic frequency components of the surface, the surface was sampled 
2500 times at At = 10 ms intervals using the CDC I 60 computer (2500 
samples is the approximate limit of the CDC 160 buffer). This block 
of 2500 was subdivided into blocks of m samples. Each block was 
Fourier analyzed using the Harm routine on the IBM 36 O computer. 

The Harm routine determines the energy of specific components of 
the signal which are given by nAf, where n = 1, 2, 3? •••• and 

Af = 1/mAt, where m must equal 2 . The highest frequency com- 
ponent that can be analyzed is given by BW = l/2At. 

Table 3.2 

Studied Sample Blocks for Harm Routine 







At 


m 


Af 




BW 




Case 


1 


10 ms 


32 


3.12 


Hz 


50 


Hz 


Case 


2 


10 ms 


64 


1.56 


Hz 


50 


Hz 


Case 


3 


10 ms 


1024 


.0975 


Hz 


50 


Hz 



50 



The surface was examined for the three cases listed in Table 



3 . 1 . The spectra obtained for the first two cases were not useful 
because they were not well enough defined over the region of in- 
terest. We are interested in the behavior of the surface between 
1 and 5 Hz. Case* 1 gave one data point in this range 3 and Case 2 
gave 3. Case 3 gave many data points between 1 and 5 Hz. 

Because Case 3 required 1024 samples , only 2 separate spectra 
could be obtained from the 2500 samples. Comparison of these two 
spectra indicated that the energy was not constant in any given 
component. Whether this energy variation was a result of second 
order interactions between the componentSj a statistical variation, 
or a combination of both was not determined. 

It was hoped that one could select a frequency component and 
analyze its amplitude variations by taking its value in each of 
the consecutive subdivisions and performing* a second order Fourier 
analysis, where AT = mAt, BW = 1/2AT. This procedure was not con- 
sidered valid for cases 1 and 2 because of the lack of definition 
in the first order spectra. Case 3 was not valid because only 2 
points were available. 

It is proposed that one could select At = 100 ms , Af == . 1 Hz, 
BW = 12.8 Hz, and m = 256 and determine 10 consecutive spectra 
from 2560 samples. This may lead to greater insight into just how 
the energy varies in and between specific components. 

The variance of the wave analyzer output was determined for 
^o ^ ^^o ~ *^35 and 5 Hz using the circuit 

shown in Figure 3-24. The 5 Hz bandwidth was obtained by replaying 
the surface analog (tape recorded # 1 ) at 3.75 ips (f^ x 10 ) and 
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Figure 3.2ii Circuit to Measure the Variance of the Output of the 
Wave Analyzer 
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setting Af^ = 50 Hz« The analyzer output was recorded on the 
second tape recorder. The gains of this recorder and the wave 
analyzer were set to one. This set the output from TR-2 equal 
to the input from TR-1, TR-2 was replayed at increased speeds 
such that the output was f^ x 1000, This output was fed to the 
Ball and Carlson sampler. The samples were recorded on the pulse 
height analyzer. These samples (Figures 3-26 through 3o28) were 
processed in the IBM 360 to obtain the variance ^ 

1 / -\2 • 2 
Var = — L (x. - X) 3 in cm , 

n 1 ^ 1 ' 

The results are plotted in Figure 3*25e 

Figure 3»25 shows the behavior of the filter. As the BW 
is decreased^ the variance approaches zero; as the BW is increased 
the variance approaches the value one would have obtained if the 
wave analyzer had not been used and the sufface had been processed 
directly. Notice, this would not be a measure of the rms height 
of the surface, but a measure of how the rms height changes. 

The PHA data, shown on Figures 3*26 through 3*28p indicate 
that the amplitude of the 3*8 Hz surface component has a stati- 
stical distribution which is close to Rayleigh for Af^ = ,03 Hz, 

As the bandwidth is increased the distribution changes from 

Rayleigh and becomes close to Gaussian for Af = ,5 Hz, 

o 

Frequency Modulation , The surface analog (f^ x 100) was fed 
into the HP 5210-A frequency meter, Figure 3.33, This meter in- 
dicates the average frequency of the signal. The needle on the 
meter varied from 3*2 to 4.8 Hz (real frequency) in an apparently 
random manner, indicating the surface is in fact frequency 
modulated , 
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Fiinction of Analyzer Band-width 




ripurr 'l.'cU PHA Data fcr Variance Measurements 
( fn = 3.6 Hz, = .03 Hz) 
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Figure 3.2? PHA Data for Variance Measurements 
( f^ = 3.8 Hz, 4fQ - .1 Hz) 
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The frequency meter has an internal discriminator network 
which gives a DC output voltage proportional to the instantaneous 
frequency* This permits measurement of frequency modulation com- 
ponents up to 25 % of the carrier frequency. For proper separation 
of the modulation from the carrier at the output, an internal low 
pass filter must be set less than, or equal to, 25% of the carrier. 
In this case a cutoff frequency of 1 Hz (real frequency) was used, 
.25 X 4 Hz. Figure 3.26 shows the input and output of the fre- 
quency meter. In this figure the lower trace is the instantaneous 
variation of the input frequency from the carrier. Figure 3*30 
shows a case where the input frequency does not vary in time. 
Figures 3.29 through 3.32 were obtained when the surface was 
filtered through the wave analyzer (f^ = 3.8 Hz, Af^ = .O3, .1, 

.5 Hz) before going to the frequency meter. 

The discriminator output, representing the frequency modu- 
lation of the surface^ was analyzed using the wave analyzer. 

Figure 3.33. The results are plotted in Figure 3.34. This is the 
spectrum of the signal which frequency modulates the surface. The 
exact mechanism which causes this modulation is unknown. The only 
thing that can be said is that it is relatively strong for modu- 
lation frequencies below .25 Hz; that is, the frequency modulation 
is very small compared to the predominate frequency in the surface. 

Conclus ions . The wind-agitated water surface is a very com- 
plicated, non-linear, dispersive medium. Results obtained from 
our analysis of the surface show that the energy in the surface 
randomly drifts through the spectrum. At one instant it is 
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.001 sec/cm 
real time 




.5 volts/cm 



Figure 3,29 Discriminator Input (upper) and Output (lower) 
for the Unfiltered Surface Analog 




,001 sec/cm 
real tim*e 



,5 volts/cm 



Figure 3.30 Discriminator Input (upper) and Output (lower) 
for Filtered Surface Analog 



( 



fo = 3.5 Hz,Zlfo = .03Hz) 
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,001 sec/cm 
real time 




.5 volts/cm 



Figure 3*31 Discriminator Input (upper) and Output (lower) 
for filtered Surface Analog 
( fo = 3.8 Hz, ^fo = .1 Hz) 




,001 sec/cm 
real time 



,3 volts/cm 



Figure 3-32 Discriminator Input (upper) and Output (lower) 
for Filtered .Surface Analog 

( f- = 3.8 Hz,4Vf„ = .5 Hz) 
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Figure 3,33 Circuit to Measure Frequency Modulation of 
the Surface 




Figure 3»3h Surface Frequency Modulation 
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concentrated at 2,8 Hz, later at 3-8 Hz, then at 3*5 Hz, and so 
on. As this energy is moving about, it is amplitude and phase 
modulating the surface components. 

Surface Wave Height and Distribution 

Electr onics . Figure 3*35 shows the circuit used to measure 
the height distribution of the wave surface. All of the components 
of this circuit have been described earlier except the sampling 
circuit and the nuclear pulse height analyzer. Both of these were 
originally used by Ball and Carlson in their measurements of the 
surface distribution. 

Calibration . To calibrate the circuit so that the channel 
numbers on the pulse height analyzer could be directly related to 
wave height, the depth of the probe was varied with the water sur- 
face smooth. For each depth the output was sampled and read on 
the pulse height analyzer. In this way the entire circuit was 
calibrated, and it was not necessary to calibrate any individual 
piece of electronic equipment. A sample calibration curve is 
given in Figure 3*36. 

Results . The wave height distribution was determined for all 
nine stations shown in Figure 3»2. The distribution at station 
1 C is given in Figure 3*37» The curves at the other stations 
were similar. The solid curve is a Gaussian curve with zero mean 
and the same variance as our data. Table 3.2 gives the values of 
rms height measured at each of the nine stations. The mean value 
of rms height was found to be 0.21 cm with a standard deviation 
of .04 cm. Because of this large variation, and because the same 
20% variation was noticed when readings were taken at one point at 
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Figure 3.36 Calibration of Surface Probe to PHA 
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Figure 3*37 Wave Height Distribution (station 1C 



different times, it was considered necessary to record the surface 
distribution and compute the rms height each time a measurement of 
the scattered pressure distribution was taken. In this way it was 
possible to ensure that the theoretical value with which we com- 
pared our data was computed using the simultaneous rms height. In 
addition to determining the variance, we measured the skewness and 
the kurtosis or **peakedness” . Skewness is defined as one-half the 
third moment divided by the third power of the standard deviation. 
And kurtosis is defined as the fourth moment divided by the fourth 
power of the standard deviation minus three divided by two. The 
skewness is less than 0.1 and the kurtosis less than 0.4 for all 
surface height distributions. For a true Gaussian distribution 
these would be zero. 



Table 3-3 



rms wave heights 



Station 



rms wave height in cm 



1 C 


0.197 


IR 


0.l6i 


IL 


0.200 


2 C 


0,228 


2 R 


0.155 


2 L 


0.263 


3 C 


0.220 


3 R 


0,183 


3 L 


0.266 
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Surface Correlation Function 



The surface spatial correlation function is defined as 



R(d) = lim C h(t,0)h(t,d)dt 



T ->oo J 

-T 



where h(t,0) is the surface height at a reference point, and 
h(t,d) is the surface height at a second point d cm from the first. 
To find this function, two probes ere placed in the surface. 
Their separation was varied using the apparatus shown in Figure 3 * 38 * 
For each separation, R(d) was computed using the circuit shown in 
Figure 3 * 39 . Oscillator frequencies of 10 and 30 kHz were used 
to reduce cross-talk between the two probes. The surface analog 
from each probe was processed in the same manner as that used to 
measure the surface spectrum. However, the DC level had to be 
more completely removed by the differential amplifier. A ,2 volt 
DC component in either signal would introduce a 4 % error in the 
correlation function. Because only one differential amplifier was 
available with an internal low pass filter, the tape recorder was 
used to filter the 10 and 30 kHz ripples from the surface analogs. 
The tape recorder operating in the FM mode has an inherent high 
quality low pass filter which can be set to 100 Hz, 1 kHz, or 
10 kHz, depending upon the speed in use. The analog computer per- 
formed the following integrations: 




T 




o 



o 



T 




o 



67 




(0 

0) 

o 

u 



68 



Figure 3*38 Variable Separation of two Surface Probes 
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Figure 3.39 Circuit to Measure Surface Correlation Function 




where A ~ auto-correlation from probe 1 
B = auto-correlation from probe 2 
C = cross-correlation between the probes 
A, B, and C were voltages proportional to the integrations. The 
absolute calibration was not important because the cross corre- 
lation function was normalized by dividing by the geometric mean 
of the auto-correlations. Integration times of 30 seconds, one 
minute, three minutes, and five minutes were tried. The normalized 
correlations function was consistent for all times. However, the 
auto-correlations displayed slow temporal and spatial variations 
of as much as 20%. It is believed that these variations are due 
to second-order interactions between waves. 

Before data were recorded at each separation, the two auto- 
correlation were simultaneously obtained. These values should 
be the same. A substantial difference indicated either a DC 
component was present in one or both of the signals, or the gains 
were different between the two circuits. These two parameters 
were adjusted until the auto-correlations were within 5 % of each 
other. An error of 5 % in the auto-correlations would be reflected 
as an error of 5% in the cross -cor relation for d less than 1 cm. 

As d is increased, C decreases because of the reduction of cross- 
correlation. However, A and B are not changed. Therefore, any 
error in A or B would cause an increased error in the normalized 
cross-correlation. This err or is variable, changing as the sepa- 
ration changes. Therefore, an error of + .1 in correlation was 
assumed for separations greater than 5 cm. 
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The correlation function was determined in the up/down-wind 
direction and the cross-wind direction. They are shown in Figures 
3.40 and 3,41. The correlation distances for these two directions 
are 3.2 cm and 10.2 cm, respectively. These results indicate that 
very little surface energy was propagated perpendicular to the 
wind direction, i.e., K » K . 

The theoretical developments consider the normalized surface 

correlation function p to go as either e or e ^ where L 

is the correlation distance. The up/down-wind function, at first 

glance, suggests a (Sin x)/x behavior. To show that it wasn’t any 

2 

of these simple relations, e ^ e ^ and [(Sin 2.2d)/3.2]/ 

(2.2d/3.2) were plotted on Figure 3*37 Tot comparisons. 

The rms slope has been shown by Beckmann to be E = ^jSr/L, 

For the up/down-wind slope this would give an rms slope of ap- 
proximately 5 ^ • 
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Figure 3»iiO Up-Dowi Wind Surface Correlation Function 
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Fig\ire Oross V.ind Surface Correlation Function 



4 • The Specularly Scattered Sound 



General 

In order to measure the effects of the surface on sound wave 
propagation, sound was directed at the surface from one transducer 
and received by another. These transducers were mounted on lo- 
cally manufactured racks which have externally adjustable angle 
dials so that the angles of incidence and reflection can be ad- 
justed to an accuracy of better than 2°, 

Spectrum 

Underwater sound at 30 kHz was specularly scattered from the 
surface described in section 3« The angles of scatter wete 0 = 21 , 
35*3 and 5l degrees. The scattered pressure was observed to be 
amplitude modulated at a low frequency, on the order of 3 Hz . 

This signal could not be directly analyzed because the narrowest 
filter on the wave analyzer was 3 Hz . At 30 kHz center frequency, 
the spectrum appeared as a pure tone; that is, the modulation fre- 
quencies were masked by the carrier. If the frequencies could be 
multiplied before analysis, thus spreading the frequencies, the 
spectrum could be measured directly, as was done in the surface 
analysis. However, the wave analyzer had a maximum frequency 
response of 50 kHz. To circumvent this problem, the envelope of 
the scattered sound pressure was detected. This demodulated signal 
was analyzed using the same techniques used for the surface ana- 
lysis, Figure 4*1* The scattered sound pressure was also analyzed 
for frequency modulation. 
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Surface 
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Ampixtude Modulation. The spectral density of the amplitude 

modulation was obtained for Q = 35 * 3 ^* This spectrum is shown in 

Figure 4 . 3 * Iri this case, the mean scattered pressure was 66 

microbars, the surface roughness was g = .l6, and the energy 

2 

density is expressed in microbars -sec. This absolute analysis 
was done for 0 = 35.3^ only. The spectra at 0 = 21° and 5l° are 
in relative units; an absolute calibration was not performed. The 
spectrum at 35 degrees was thoroughly analyzed to find its re- 
lationship to the surface spectrum; the spectra at 21 and 51 degrees 
were obtained to confirm the findings. 

Frequency Modula t ion . The circuit shown in Figure 4*5 was 
used to measure the frequency modulation of the scattered sound 
pressure. To ensure accuracy of the experiment, a frequency syn- 
thesizer was used at a frequency of 30j000*0 *1 Hz. To analyze 

the received signal, the HP 5210 frequency meter, which was des- 
cribed in section 3? was used. The cutoff frequency of the output 
filter (1 kHz in this case) must be chosen such that the carrier 
does not introduce ripple in the output. When a pure tone is 
put into the meter, the output is a DC voltage between 0 and 1 volt. 
The actual level depends on the frequency and meter range in use. 
When the carrier is frequency modulated, the output varies about 
this reference DC. The amount of variation is proportional to the 
frequency deviation, and the frequency of the variation is equal 
to the modulating frequency. 

In order to calibrate the system, sound pressure was specularly 
scattered from a smooth surface (0 = 35 . 3 °). The synthesizer fre^ 
quency was changed ^ 10 Hz and the output DC voltage was found to 
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Figure U.2 Sound Pressure Amplitude Modulation Energy Spectral Density, G =2l' 
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Figure ii.i^ Sound Pressure Amplitude Modulation Energy Spectral Density, 
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deviate + 10 mv. This procedure was repeated for several frequency 
changes down to . 5 Hz . When the frequency change was 1 Hz , 
a deviation of 1 mv was detected in the output. However, for 
frequency changes below 1 Hz no deviation was detectable above the 
noise. These deviations* were in agreement with values calculated 
following the instructions in the HP 5210 manual. These results 
indicated the frequency deviation in the analyzed signal must be 
greater than 1 Hz to be detected. 

The surface was roughened and the sound scattered from it. 

It was impossible to detect any change in the output from the fre- 
quency meter. This result led to the following conclusions. 1) If 
there was a frequency shift in the carrier it was less than 1 Hz. 

2) If there was any frequency modulation the magnitude of the de- 
viations was less than 1 Hz. 

Analysis . Examination of Figure 4.3 leads to the following. 

At Q = 35.3 degrees the spectrum of the amplitude modulation had 
a form which closely resembled the surface spectrum, except for 
the dips near the peak and the added peak at about 6 Hz. The 
problem was to find the reasons for this behavior. The first 
thought was doppler shifting of the pressure upon scattering, al- 
though Parkins and Wetzel had shown that doppler shifting should 
not occur in the specular direction. 

It was discovered that this spectrum was in fact a result of 
doppler shifting. Figure 4.6 depicts the scattering geometry in 
the vertical plane of incidence defined by the axes of the two trans- 
ducers and the perpendicular to the surface at the center of the 
scattering area. Consider an object as shown. When sound pressure 
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Figure Ii,6 Scattering Geometry 



is incident on this object, the majority of the energy will be 
scattered along path B. However, there will be diffraction side 
lobes to this scattered pressure. 

Consider the case in which the object moves with velocity V 
and, at a given instant is positioned such that the side lobe is 
coincident with the receiving transducer. The pressure in this 
lobe can be shown to be doppler shifted by the scattering object. 
From geometry the amount of doppler shift is given by 

^ [Sin(0 - p)- sin (0 + o')] 

Considering O' — - p and both small, this equation reduces to 

^ (P - Q')Cos 0 

f d - r Cos 0 

Note, for specular scattering, p = 0 and f ^ = 0; that is, this 
agrees with Parkins and Wetzel. 

Wetzel and Parkins, considering a resonant interaction between 
the acoustic wave length \ and the surface wave length A, show the 
diffraction grating equation is 

p = ^ (Sin 0^ - Sin 0^) 

where p is the order of the lobe, 0^ is the angle of incidence and 
0^ is the angle of reflection. Letting = Q - p , and 0^ = Q + O', 
and considering p = 1; that is the first side lobe, this equation 

reduces to . 

A = 

2p Cos 0 
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If V = considered the velocity of the diffraction grating, 

that is, the velocity of the surface wave component having a wave 
length A, 

= f , r 2B Cos 0 
d surf X ^ . 

which becomes f , = f ^ when A = X/26 Cos 0 is inserted, 
d surf ^ 

This development shows that the dopplex frequency shift is 
equal to the frequency of the surface component causing the shift. 

To have a shifted component present in the received sound certain 
conditions must be met. There must be a path BB’ which satisfies 
the diffraction equation in such a way that the first side lobe 
falls on the receiver. This implies that for each doppler shifted 
component there is an appropriate corresponding angle p in the plane 
of incidence (and similar angles outside this plane) . No more than 
one shifted component can arrive along a given BB’ path. To change 
fd> ^suixf change; then the relationship A = X/2p Cos 0 no longer 

holds; and the side lobe moves off the receiver. 

When two tones arrive at the transducer slightly different in 
frequency, they combine to produce amplitude modulation at the 
beat frequency. The frequency of the combined signal is somewhere 
between the two tones, depending on their relative strength. The 
model being considered deals with the spectrum of th^ amplitude 
modulation, that is, the spectrum of the beat frequencies. A given 
beat frequency equals the amount of doppler shift f^. 

As reported, it was not possible to determine directly the 
frequency spreading of the carrier which resulted from this inter- 
action; the measurement of the amplitude modulation spectrum was 
conducted in place of the direct measurement. 
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The dips near the peak of the spectrum in Figure 4*8 sire 
difficult to explain. They seem to be a result of interference. 

For the geometry used in this experiment, these dips occur when 
there is a (n + j n = 1, 2, ••• path length difference between 

paths and AA. For f^ = 2.6 Hz this difference is ; for 

^d ~ difference is 3/2X (X = 5 cm). This relationship 

would explain destructive interference if the signals arriving 
along the two paths were coherent. However, it has been shown 
that they are incoherent. There is a possibility that there are 
surface elements outside of the plane of incidence which satisfy 
both the angle relation (for the proper frequency) and the path 
difference requirement for interference. The geometrical relations 
necessary to test this speculation have not been developed. 

Another possible explanation for the peaks and troughs is that 
they are due to higher diffraction orders. 

The secondary peak that occurs at f^ = 6 Hz, Figure 4«8 is 

attributed to the effect of the transducer side lobe. This is 

discussed under results and conclusions. 

Study of 

Results and Conclusions . A The surface energy spectral density 
E(f^), the scattered sound amplitude modulation energy spectral 
density £(f^), the sound pressure level at the surface where the 
path BB’ scatters, SPL(p), shows that the shape of 
resembles the shape of E(f^), but is displaced; the peak of ^(f^) 
moves to the left (lower frequency) as 0 increases, as predicted 
by the doppler shift equation f ^ ^ 2B Cos 0; for a given fre- 

quency, the SPL(p) becomes smaller as 0 increases, again as pre- 
dicted by the doppler shift equation; and the secondary peaks in 
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Frequency (Hz) 

Figure U.8 Relationships of E(f^), SPL(^) and yS , for 0=35.3 



G(f^) at 6 Hz, 0 = 21 and 35.3 degrees, correspond to doppler 
shifted energy arriving at an angle which corresponds to the side 
lobe of the transducer radiation pattern. Not shown on these 
figures, but found from the sound pressure magnitude studies, is 
the fact that as 0 increases up to about 75 degrees, the variance 
of the amplitude modulation not only moves to lower frequencies, 
but also decreases in magnitude as Q increases. 

An explanation of how the doppler shifted energy is produced 
and arrives at the transducer was presented in the analysis section 
however, no attempt was made to explain the variation of the mag- 
nitude of the energy. It is contended that the energy in a given 
doppler shifted component is a function of both the energy con- 
tained in the surface component that causes the shift and the 
sound pressure level incident on the surface at the point where 
the shift originates. 

fi(f^) ^ F E(f^) 

where F = F[SPL(p) ,A] 

A is the surface area that satisfies the diffraction condition. 

It must be emphasized that the explanation of the amplitude 
modulation spectrum has been based on the simplifying assumptions 
that (a) scattering comes principally from the water surface close 
to the plane of incidence and (b) the diffraction grating equation 
can be applied to the very limited surface area where the rays are 
incident at angle p ^ A where A « p defines a bundle of rays of 
common angle of incidence. 
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After completion of this research, a recent paper by Marsh 
and Mellen (I 966 ) was discovered in which the amplitude modulation 
in specularly scattered sound was ascribed (without proof or 
detailed data) to interaction between the unshifted and the doppler 
shifted energy. 

Distribution of Scattered Sound Pressure 

Underwater sound at frequencies between 20 kHz and 100 kHz 
was specularly scattered from the model sea surface described in 
section 3 at angles of incidence between 35 and 85 degrees, with 
particular emphasis on angles greater than 70 degrees. 

The circuit diagram for measuring the distribution of the 
scattered sound pressure is shown in Figure 4.10. This circuit 
yields a distribution from the pulse height analyzer. When the 
calibration is applied, the distribution of sound pressure is 
obtained. 

Calibration . The calibration is done in still water. In 
order to obtain the proper calibration, the output voltage of the 
receiving transducer is amplified a fixed amount (usually X2 or 
X5) by an HP 46?A amplifier. The output of this amplifier is 
then measured using a vacuum tube voltmeter. The signal is sampled 
and measured by the pulse height analyzer . This determines the 
mirror reflection value of mean square pressure. Then the selector 
switch on the HP 46?A is changed to variable and the output changed, 
read on the voltmeter, and read on the PHA. This calibrates the 
PHA to the output of the amplifier . The amplifier is returned to 
its original fixed amplification setting, and the calibration is 
complete. The resulting calibration line is a plot of transducer 
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Figure U.IO Circuit for Measuring Distribution of Scattered Sound Pressure 




output times a multiplication factor (two or five) versus channel 
number on the PHA. A sample calibration curve is given in Figure 
4.11. To convert this calibration from output voltage to sound 
pressure, the response curve for the transducer was used* This 
curve gives the output in db re 1 volt/microbar as a function of 
frequency. 

Measurements . With the surface roughened as described in 
section 3, the received sound signal was sampled at a sampling 
rate of one thousand samples per second, and the samples were 
recorded on the PHA for seven minutes. In all, seventy-four sets 
of data were recorded and are provided in Table 4.1. For g less 
than 0*3 the theoretical value of relative mean square pressure 
is e ^ and it was this theoretical value with which we compared 
our data. Except for four sets of data, all of our g^s were less 
than 0.23. 

As mentioned in section 3? during the same time that sound 
pressure data were taken, the wave height analog was recorded. 

This signal was processed later to determine the rms wave height 
which existed during the period of the sound scattering experiment. 

A program (Appendix A) was used in the IBM 36O computer to 
provide, given the PHA distributions for the wave surface analog 
and the sound pressure analog and the necessary calibration factors, 
the proper values of the rms wave height, the skewness and kurtosis 
of the wave height distribution, the mean scattered pressure, the 
variance, skewness, and kurtosis of the sound pressure distribution, 
the roughness factor g, the theoretical value of relative mean 
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Table 4. 1 Scattered Sound Data 
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Table 4.1 (Continued) 
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square pressure e and the ratio of measured relative mean 
square pressure to e Figure 4-12 shows an example of the 

distribution of the specularly scattered sound pressure. 

Results . Figure 4.13 is a graph showing the ratio of the 
measured relative mean square pressure to the theoretical value 
e ^ as a function of the angle of incidence for various frequencies 
for g < 0.2. The dashed line at 1 corresponds to the theoretical 
value e The solid line is the theoretical value including the 

shadowing function of Wagner. 

Figure 4*l4 shows the variance of the sound pressure as a 
function of the cosine of the angle of incidence for various 
frequencies . 

Figure 4.l6 shows the variance at a particular angle as a 
function of frequency. 

Conclus ions . Figure 4*13 clearly shows that for the largest 
angles of incidence (smallest grazing angles) the measured values 
of relative mean square pressure are lower than the value predicted 
by the theory of Beckmann and higher than predicted using the 
shadowing function of Wagner. It seems clear that shadowing does 
occur; however, the calculation of Wagner over corrects for the 
shadowing by ignoring multiple scattering and the wave nature of 
sound. Wagner’s theory uses ray geometry and assumes that any 
particular point on the surface is either illuminated or not 
illuminated. Because the wavelength of the sound is the same 
order of magnitude as the rms height of the surface, it should 
be expected that bending in the wave front will occur and areas 
which Wagner predicts will not be illuminated will in fact be 
partially illuminated . 



102 




103 



Figure U.12 Distribution of Specularly Scattered Sound Pressure (© = 77.6°, f= 85 kHz) 
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Variance 




Figiire U.l5 Variance vs. Angle of Incidence for Various 
Frequencies 
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Variance 




Frequency (kHz) 

Figure U.16 Variance vs. Frequency for 72° Angle of 
Incidence 
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The measured variances were all equal to or greater than the 
variances predicted by Beckmann. This is to be expected. If 
independent random variables are addedj, the variance of the sum 
is the sum of the variances. Therefore, it is believed that the 
variance predicted by Beckmann is the minimum variance to be 
expected with any additional variances due to the electronics, 
shadowing, and multiple scattering adding to this predicted value. 
The dependence of the variance on cos^ 0 is seen in Figure 4*l4, 
for angles less than 75 degrees. This dependence is in agreement 
with Beckmann and indicates that variances introduced by the ex- 
periment are nearly constant. 

Also it can be seen in the same figure and also in Figure 4.15 
that the variance increases rapidly for angles greater than 75 
degrees. It is believed that an additional variance due to the 
shadowing causes this behavior. As shown in section 3, the 
Gaussian distributed rms slope is approximately 5 degrees. This 
implies that 99% of the slopes will be less than three times this 
value or 15^. Therefore, when the angle of incidence is approxi- 
mately 75^ (grazing angle of 15^) p some of the surface should begin 
to be shadowed. The variance in the shadowing effect is then added 
to the other variances. 

Figure 4.l6 shows the variation as a function of frequency 
for a particular angle of incidence (72^). It can be seen from 
this graph that the variance is approximately proportional to 
the second power of frequency. This is not in agreement with the 
prediction of Beckmann that the proportionality should be with the 
sixth power of frequency. No explanation for this discrepancy has 
been found. 



108 



It is interesting to notice that the various theories of 
sound scattering always assume a Gaussian distribution of surface 
heights and, for g « 1, predict a Gaussian distribution for the 
scattered sound. In fact, neither the real sea nor our model 
shows an ideal Gaussian distribution; and, of course, the scat- 
tered pressure is not an ideal Gaussian either. A virtue of 
the model experiment is that we can measure the resulting sound 
scatter distribution for a height distribution that shows re- 
alistic non-zero skewness and peakedness. 
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5 . Summai i 

Experiments were conducted to measure the distribution and 
spectrum of sound which was specularly scattered from a wind gen- 
erated model sea surface* The energy spectral density of the 

"9 • 3 

surface peaks at 3*8 Hz, dropping off at f * with increased 
frequency near the peak, becoming f for frequencies greater 

than 5 Hz. This frequency behavior is in agreement with the 
Neumann-Pier son spectrum and Phillips* theory which predicts f ^ 
dependence. The rms height of the surface was about 0.21 cm. 

The wave height distribution was nearly Gaussian. The spatial 
correlation function of the surface was determined for both up/ 
down-wind and cross-wind directions. The correlation distances 
were 3.2 and 10.2 cm, respectively. From the correlation dis- 
tance and the rms height, the rms slopes were calculated for up/ 
down* and cross-wind directions. These values were .OO 96 and 
. 0008 , respectively. These results confirmed the assumption that 
the surface propagation vectors were nearly parallel down-wind 
with very little divergence. 

The distribution of the specularly scattered sound pressure 
was measured for various angles of incidence from 35 ^ to grazing 
and various frequencies between 20 and lOO kHz, with particular 
emphasis on angles near grazing. Comparison was made with theories 
that predict the values of the mean and variance of the specularly 
scattered sound distribution. The mean values are in agreement 
with theory except at very near grazing indidence where shadowing 
causes less specularly scattered energy than predicted for no 
shadowing and more than predicted for geometrical shadowing. The 
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measured variances were generally greater than predicted due to 
the additative experimental variances. The observed dependence 
of variance on the fourth power of the cosine of the angle of in- 
cidence was in good agreement with theory. However, the variance 
was proportional to the second power of frequency rather than the 
sixth power as predicted. 

The spectrum of the specular scattered sound contains a 
predominant component at the frequency of the radiated sound. 

In addition it contains doppler shifted components, although 
theory predicts that the spectrum of purely specularly scattered 
sound should not contain doppler shifted components. 

The doppler shifted energy arrives at the leceiver from 
directions other than the purely specular. This is a result 
of finite divergence of the incident sound. It was found that 
the doppler shifted spectrum is directly related to the frequency 
spectrum of the surface. The energy level of the doppler shifted 
components is a function of the surface energy spectral density 
and the incident sound pressure at the off-axis points. 
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